Role of bound pairs in the optical properties of highly excited semiconductors: a self 

consistent ladder approximation approach. 
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Presence of bound pairs (excitons) in a low-temperature electron-hole plasma is accounted for by 
including correlation between fermions at the ladder level. Using a simplified one-dimensional model 
with on-site Coulomb interaction, we calculate the one-particle self-energies, chemical potential, and 
optical response. The results are compared to those obtained in the Born approximation, which 
does not account for bound pairs. In the self-consistent ladder approximation the self-energy and 
spectral function show a characteristic correlation peak at the exciton energy for low temperature 
and density. In this regime the Born approximation overestimates the chemical potential. Provided 
the appropriate vertex correction in the interaction with the photon is included, both ladder and 
Born approximations reproduce the excitonic and free pair optical absorption at low density, and the 
disappearance of the exciton absorption peak at larger density. However, lineshapes and energy shifts 
with density of the absorption and photoluminescence peaks are drastically different. In particular, 
the photoluminescence emission peak is much more stable in the ladder approximation. At low 
temperature and density a sizeable optical gain is produced in both approximations just below the 
excitonic peak, however this gain shows unphysical features in the Born approximation. We conclude 
that at low density and temperature it is fundamental to take into account the existence of bound 
pairs in the electron-hole plasma for the calculation of its optical and thermodynamic properties. 
Other approximations that fail to do so are intrinsically unphysical in this regime, and for example 
are not suitable to address the problem of excitonic lasing. 



I. INTRODUCTION 

The binding of a gas of oppositely charged fermions 
into a bosonic gas of bound pairs has a fundamental in- 
terest in condensed matter physics, in problems ranging 
from the state of hydrogen to the nature of the electron- 
hole plasma in Si and Ge. This pairing takes its ori- 
gin from the Coulomb correlation between the charged 
fermions, and strongly affects the statistical and thermo- 
dynamic properties of these systems. In semiconductor 
physics, the role of the Coulomb correlation and pairing 
in the description of the electron-hole plasma in Si and 
Gc has been extensively investigated. Refined descrip- 
tions of both the ground state and the thermodynamic 
properties of this system have been developed in the past 
thirty years. In particular, an excitonic insulator ground 
state was first proposed and then shown to be less stable 
than a simpler electrcm-hole liquid in Si and Ge, due to 
band-structure effects^ Thermodynamics of the electroo.- 
hole gas in sinroler crystals was also addressed by Hauga, 
and Zhu et alB aiming at deriving a complete phase di- 
agram from variational approaches. 

In this paper we will deal with direct gap semicon- 
ductors and we will discuss the thermodynamics of the 
photoexcited electron and holes at quasi-thermal equilib- 
rium. Our main concern is to determine how the inclu- 
sion of pairing effects in the theoretical description, as 
done for the systems mentioned above, can affect pre- 



dictions on the shape of optical spectra. In the optical 
properties of semiconductors, the effect of the Coulomb 
electron-hole correlation comes into play to explain the 
characteristic excitonic absorption and emission. These 
excitonic features can be described by including vertex 
corrections in thiSj interaction of the electron-hole pair 
with the photon.cl In particular, excitonic absorption is 
found in an empty crystal, when a background electron- 
hole plasma is absent. However, when carriers are electri- 
cally or optically injected, they reach quasi-thermal equi- 
librium through scattering, and a background electron- 
hole plasma is eventually formed. In this typical situa- 
tion the absorption is modified, and spontaneous emis- 
sion (photoluminescence) also takes place. Absorption 
and emission give valuable information about the state 
of the plasma, but a theory describing both the ther- 
modynamics and optical properties of the plasma is re- 
quired in order to extract this information. Semicon- 
ductor Bloch equations, and their evolutions, are cer- 
tainly among the most known and used theories for this 
purpose.cl In these theories, the electron-hole plasraa has 
been originally treated at the Hartree-Fock leveljj and 
eventually screening of the Coulomb interaction was also 
included at various levels of approximation. Q However, 
none of these refinements accounts for the existence of 
bound excitons in the plasma. This is certainly a weak 
point at low temperatures, when simple thermodynamic 
arguments show that condensation of electron-hole gas 



into a gas of bound-electron hole pairs is favorable in a 
wide range of densities. 

In the following sections we will show the relevant 
changes in the optical properties of the semiconductor 
when the existence of excitons at low temperatures is 
taken into account. This is clearly a leap beyond the 
standard description of the interaction of light with the 
semiconductor as described above. Our purpose is quali- 
tative, but we can still gain solid understanding of the un- 
derlying physics. In our opinion at this stage of develop- 
ment it is not wise to aim at improving also quantitatively 
existing theories which are already far too involved. We 
keep the problem simple enough so as to always control 
and understand the resulting physics. We thus consider 
an accurately solvable model, where the system is one- 
dimensional. Coulomb interaction is on-site, and electron 
and holes are spin-polarized. Excitonic correlation is also 
expected to be stronger in one dimension than higher di- 
mensions. Moreover, for the single-pair problem, exactly 
a single bound state is obtained with the on-site interac- 
tion. In our description of the electron-hole plasma, we 
either include carrier scattering at the Born level, which 
does not describe bound pairs in the plasma, or we in- 
clude enough correlation so as to describe the pairing 
into bound states. For this purpose—ladder diagrams in 
the electron-hole scattering kernel3~B are included, aud, 
the electron and hole self energies calculated with it.t2l 
The resulting Dyson equations for the Green functions 
are solved self-consistently. The approximation is there- 
fore called self-consistent ladder approximation (SCLA). 
Screening is expected to be weak in one dimension,0 due 
to reduced screening phase-space, and we neglect it al- 
together. We remark that in the SCLA we calculate the 
single particle propagators and self-energies only, thus in 
this sense, the theory remains purely fermionic. We do 
not map the model onto bound and un-bound states at 
any time. However, it is possible to show analytically 
that at low temperatures and densities, a bosonic gas of 
excitons is effectively described, and that self-consistency 
effectively introduces scattering between these excitons. 
When this analogy is extended to larger densities and 
temperatures, where bound and free carriers are expected 
to coexist, we understand that self-consistency also ac- 
counts for multiple exciton-free carrier, and free carrier- 
free carrier scattering. Thus, a large amount of correla- 
tion, well beyond that described in simpler Born approx- 
imations is included in the SCLA. We remark that we 
also assume quasi-thermal equilibrium in the electron- 
hole gas, which is well justified as the characteristic ra- 
diative recombination (several hundreds of picoseconds) 
is much slower than the typical scattering time in the 
density range considered. The theory is however written 
so as to be readily extended to the non-equilibrium situ- 
ation. We finally show how to correctly calculate photon 
absorption and emission (photoluminescence or PL) in a 
conserving sense, i.e. respecting f-sum rules. We only 
neglect pol^aritonic effects, which are very weak in one 
dimension. t3 



Within the considered models, we highlight important 
trends in optical and thermodynamical properties. First, 
Born approximations overestimate the chemical poten- 
tial at low temperatures when a consistent fraction of 
pairs is bound into excitons. Second, energetic stability 
of the excitonic emission peak as a function of density 
is better described in the SCLA, even well beyond dis- 
appearance of excitonic absorption peak at large excita- 
tion densities. Third, the large excitonic gain at densities 
just below those where excitonic absorption disappears is 
clearly incorrectly described in simpler Born approxima- 
tions. As these trends have a clear physical origin, they 
are expected to hold qualitatively even when more realis- 
tic interaction potentials (eventually including screening) 
are considered. We conclude that these results do have 
important experimental and theoretical relevance, signal- 
ing the necessity of an adequate description of the excited 
semiconductor in the low-temperature region. For exam- 
ple, current description of the electron-hole plasma in 
the semiconductor Bloch equations is clearly insufficient 
to address excitonic (inversion-less) lasing. 

The paper is organized as follows. In Sec. O we intro- 
duce and explain how to implement the self consistent 
ladder approximation, and the Markov-Born approxima- 
tion and self-consistent Born approximation. Single par- 



[II, where we show 



tide properties are analyzed in Sec. 
how excitonic correlation appears in the electron and hole 
spectral functions. Chemical potentials and self-energies 
are then compared. In Sec. IV we analytically show how 



the SCLA at low temperature and density describes an 
interacting exciton gas. Optical properties of a semicon- 
ductor quantum wire are calculated in Sec. |V| for the 
three approximations. Absorption and emission as func- 
tions of the density are also presented, and discussed in 
this section. Conclusions are drawn is Sec. 
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II. BORN AND SELF CONSISTENT LADDER 
APPROXIMATION 

In this section we introduce the SCLA, the simpler 
Markov-Born approximation (MBA) and self-consistent 
Born approximation (SCBA). We assume thermal equi- 
librium in the electron- hole plasma, but, in order to avoid 
analytic continuation problems, we dcujiot resort to the 
Matsubara Green function technique jl3 which becomes 
numerically delicate when the spectral function shows a 
large number of poles. Instead, we work directly in the 
real time/frequency space. This has the additional ad- 
vantage that the_approach can be readily extended to 
non-equilibrium.EJpWe thus consider four single particle 
Green's Functions:t3 the retarded, advanced, lesser and 
greater Green's functions defined as 

G+(l, 2) = -leih ~ t2)(T[c(l)ct(2)]+), (la) 

G-(l,2)-j0(t2-ti)(r[c(l)ct(2)]+>, (lb) 

G<(l,2)=z(ct(2)c(l)), (Ic) 



G>(l,2) = -^(c(l)ct(2)). 



(Id) 



Here T is the time ordering operator, 1 = (xijii, cri), 
where xi is the position, ii the time, and cti is the spin 
index. c(l) is the electron annihilation operator and [ ]+ 
is the anticomniutator. Similar expressions for the hole 
Green's functions hold, with d{l) and d){l) the hole an- 
nihilation and creation respectively. Only two of the 
four functions above are independent, and in particu- 
lar, we will consider the retarded and lesser (or corre- 
lation) functions. Formally, a single Green's function 
G(xi, Ti, 2:2, T2), may be introduced for compactness of 
notation, with ti and T2 defined on the Keldysh contour, 
having the two branches shown in Fig. ^. Following sim- 
ple rules, any equation with times defined on the Keldysh 
contour can be written back in a set of equations defined 
on the ordinary time axis for the above four Green func- 
tions (see for instance Rcf. [l5| ). 
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FIG. 1. The Keldysh contour. Upper branch, normal time 
ordering, lower branch, reversed time ordering. 

The statistical averages are done in the macro- 
canonical ensemble: (...) — Tr{p...} with p — 

1/T and H is the total Hamiltonian of the interact- 
ing electron-hole system. In principle, electrons and 
holes have independent chemical potentials, and densi- 
ties. However, we are interested in the description of 
a laser excited semiconductor where the density of elec- 
trons and holes is the same, and charge is balanced. We 
will also assume that the temperature of the two com- 
ponents of the interacting gas is the same. In fact, after 
excitation, thcrmalization of carriers is mainly driven by 
exchange of energy with the phonon thermal bath which, 
at equilibrium, leads to T^ = Th = T . Moreover, for sim- 
plicity we consider the same electron and hole masses, 
giving fj,e = Hh = ^J■ and the same Green functions and 
self-energies. Extension of the theory to different masses 
is straightforward. Finally, the Green functions above 
depend only on the time difference t2 — ti as the system 
is stationary, and on the relative distance X2 — xi as the 
system is homogeneous. 

In this thermal regime, the retarded and lesser Green 
functions defined above are, also related by the Kubo 
Martin Schwinger relationslij 

G<(k,uj) = -2^[G+{k,uj)]f{uj-n) (2a) 

G>{k,Lu) - 2Q[G+(fc,u;)](l - f{uj - /i)), (2b) 

where we considered the Fourier transforms with respect 
to the relative time ti— ^2, and with respect to the relative 
position a;i—a:;2. The function /(a;) = [ea;p(cij/T) + l]^^ is 



the Fermi function. In the stationary case G~^{k,u},a) = 
G~ {k,uj,a) contain the information on the spectral 
properties of the quasiparticles, given by the one parti- 
cle spectral function A{k, oj, a) = — 2S5[G+(fc, lj, a)]. The 
correlation functions instead contain information on the 
quasi particle occupation number G^{k, ti — t2 ^ t, a) = 
N{k,t,a), and G>(fc,ti = t2 = t,cr) = 1 - N{k,t,a), 
where N(k, ta) are the occupation numbers. 
The on-site Coulomb interaction reads: 



V{x) = ±a5{x) , 



(3) 



the potential is repulsive for electron-electron and hole- 
hole, and attractive for electron-hole interaction. An at- 
tractive delta-like potential between the electrons and the 
holes gives one bound state only for the pair problem, 
and a can be chosen such to reproduce the exciton bind- 
ing energy measured in the limit of low density. The 
constant a has the dimensions of energy x length, the 
exciton Bohr radius is given by a^ = ?i /{ma) and the 
exciton binding energy is Ei, — ma? /2h where m is the 
reduced mass of the electron-hole system. We will use 
the units U — a — m — 1 throughout the paper. In these 
units as = 1 and Ei, ~ 0.5. We remind that typical val- 
ues for semiconductor quantum wires are qb = 10~^ cm 
and Et = 10 meV.EII With the above effective potential 
of Eq. ^ the interaction Hamiltonian is 



^c-^E^^(l)^^(l')c(l')c(l) 

1 

-5]ct(l)dt(l')c(l)d(l'). 



(4) 



Here 1' = {xi,ti,a' ^ a) for the electron-electron 
and hole hole interaction, while 1' = {xi,ti,a') for the 
electron-hole interaction, as due to the fermionic nature 
of the carriers, only interaction with opposite spin is al- 
lowed in the intraband term for a contact potential. In 
the following, we assume a spin-polarized system of elec- 
trons and holes, so that we keep only the last term de- 
scribing the electron-hole interaction in the Hamiltonian 
(Q). Generalization to both spins is straightforward. 

The Dyson equation for the single particle G(l,2) 
propagators read 

G(l, 2) = Go(l, 2)+ fdS d4Go(l, 3)I](3, 4)G(4, 2) , (5) 

where time integrations are over the Keldysh contour. 
Go is the free propagator. The functional dependence of 
the self energy E on the single particle Green's functions 
determines the degree of approximation. In particular, 
it can be expressed through an electron-hole scattering 
kernelT(l,2;l',2'): 



S](l, 2)^i f dl' d2' T{1, 1'; 2, 2')G{2\ l') 



(6) 



In the Born approximation, 

r(i,i';2,2') = y(i,i')+ 

+ V{l,l')iG{l,2)G{r,2')V{2,2'), (7) 

where 1^(1, 1') = V{xi — xi')6{ti, ty), and 6{t, t') is the 
Dirac deha function extended on the Keldysh contour. 
The self-energy I]{k,uj) can be calculated at the single 
particle pole uj = k'^/A. In this case we are neglecting 
memory terms in the scattering and thus using a Markov 
approximation; this corresponds to the MBA. In the 
SCBA, the same scattering kernel of Eq. ^ is used, but 
the full energy dependent structure of the resulting self- 
energy given in Eq. (^) is used in the Dyson equations 
Eq. (ph. The problem is then solved self-consistently, 
as explained later. In this case, memory effects, be- 
yond the Markov approximation, are also included. In 
the SCLA, the kernel r(l, 1'; 2, 2') is the solution of the 
Bethe Salpeter equations in the ladder approximation, 
which for a generic potential reads 

r(l,l';2,2') = y(l,l') + 



d3 d'i' y(l, l')JG(l, 3)G(1', 3')r(3, 3'; 2, 2') (8) 



For a 5-like potential, V{\, 1') = —(5(1, 1') and T depends 
only on 2 — 1. The reduced equation reads: 



(9) 



T(l; 2) = -l-i d3 G{1,3)G{1, 3)T(3; 2) 



The Eqs. (H), (g) and (||) are schematically represented 
in Fig. I for both the SCBA and the SCLA. 
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FIG. 2. (a) The Dyson equation, thick lines are the dressed 
one-particle Green function, thin lines the bare ones, (b) The 
self-energy with the scattering kernel T{1;2), for the on-site 
Coulomb potential, (c) The scattering kernel T(l;2) for the 
Born approximations, (d) The scattering kernel r(l; 2) in the 
ladder approximation, given by the Bethe-Salpeter equation. 



In the self-consistent approximation we solve self- 
consistently Eqs.(^ and (|g) with the scattering kernels 
given by Eq. (^ or (Q) in the Born or SCLA case respec- 
tively. We explain in details how numerical self consis- 
tency for a fixed temperature and chemical potential /i 
has been implemented. 

(i) We start with G^{k,uj) of the form G+o — (^ — 
fc2/4-hi7)-\andG<o(fc,w) = -25[G+o(fc,w)] /(w-^)- 

(ii) In order to take advantage of the local character 
of the electron-hole interaction, we use the Fourier trans- 
form of the Green functions of Eq. (^, G(l,2), which 
actually have only a spatial dependence on r = xi — X2 
as the system is homogeneous. We calculate the free pair 
propagator defined as 



i7(l;2) = *G(l,2)G(l,2) 



(10) 



The retarded and correlation parts are calculated passing 
from the Keldysh contours to the usual time axis: 



H<{r,t) ^ iG<{r,t)G<{r,t), 



(11) 



H+{r, t) = iG+{r, t)G+{r, t) + 2iG+(r, t)G< (r, t). (12) 

The dependence on the relative distance has been explic- 
itly shown. The dependence on the relative propagation 
time t — t2 — ti stems from the stationary condition. 
The pair propagator is then Fourier transformed to (fc, w) 
space. 

(iii.a) For the Born approximations, we use Eq. (M) 



T<{q,u;)^H<{q,u;) 



(13a) 
(13b) 



(iii.b) For the SCLA case, Eq. (||) is readily solved as: 



T+{q,Uj) = -{l + H+iq,Lu))-\ 
T<{q,u;)^\T+{q,cj)\'H<{q,cj). 



(14a) 
(14b) 



(iv) The electron self-energy is then calculated in real 
space as in Eq. m). The retarded and correlation func- 
tions, defined on the usual time axis, read: 

I]<(r,i) = iT<{r,t)G>{-r,-t) = -iT<{r,t)[G> {r,t)]* , 



S+(r, t) = -iT<{r, t)G^{-r, -t)+iT+{r, t)G< (-r, -t) 



= -tT<{r,t)[G+ir,t)r+iT+ir,t)[G<ir,t)r. 

(v) The new electron Green's functions are then calcu- 
lated as 

G+(fc, uj)^{uj- k^/A + S+(fc, oj) + i7)-\ 



G<(k,u)^-2^[G+{k,uj)] f{uj~n). 



The procedure is repeated through step (ii) until self- 
consistence is reached. Fast Fourier transforms are per- 
formed on a finite grid of 16384x128 points for the fre- 
quency and wavevector domains respectively. The fre- 
quency and wavevector domain (fc,w) is (-20,20)X(-3,3). 
Due to the finite k-range, we obtain Ei, — 0.4. The exter- 
nal 7 > in (v) is adiabatically switched off during self 
consistency. In this way the imaginary part of the final 
G is provided by the interaction only. For the non-self- 
consistent Born approximation (or Markov-Born) , we use 
the frequency independent, on-pole S^((jj — fc^/4) in the 
Green function, and stop the procedure at step (v). 

All of the considered approximations are conserving 
in the sense of Kadanoff and Baymila the total density, 
total momentum, total energy and total angular momen- 
tum of the system are conserved. This is a fundamental 
property of any approximation for the self-energy, other- 
wise unphysical results may result. We will come back to 
this point again in the calculation of the optical response 
of the system in the various approximations. 

The polarized electron-hole gas with a contact poten- 
tial interaction that we consider here can be mapped 
onto a single band Hubbard model with spin 1/2 and 
attractive interaction, in the limit of infinite width of the 
band, and infinite on-site interaction so as to produce fi- 
nite Eb (continuum limit). For the Hubbard model with 
repulsive interaction, the self-consistent ladder approxi- 
mation is-known as fluctuation exchange approximation 
(FLEX)JIj Even though in this repulsive case the opening 
of an Hubbard gap is not reproduced, in the attractive 
case the gap is of different nature andis notoriously well 
described within this approximation.E2l 



III. SINGLE PARTICLE PROPERTIES 
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FIG. 3. Electron spectral function at fc = as a function 
of the total carrier density. T — 0.1. 















n=0.001 








n=0.004 


- 








n=0.05 




.J 




n=0.1 








: 




-' .^ ••■■ / ^^y \\ ■'■■■. ~^-^~-.^ ; 


■^'^ ■■■■■ /^A/ V--^^'''----- ^"^- 










^"Z^^""""'^ ^"""--'^ ^^ 






^^^^ 1 . 1 . 1 . 1 . 1 ^^^ 



We plot in Fig. || the electron spectral function at 
fc = for different densities, obtained in the SCLA. The 
temperature T = 0.1 ^ £'b = 0.5. For this particular fig- 
ure we used the k-range (-6,6) in the numerical solution. 
For n ~ 0.001 the electron spectral function has a very 
narrow Lorentzian shape. For larger densities, a satellite 
structure appears in the low energy side of the spectral 
function. This structure is located at the exciton bind- 
ing energy, below the main peak, and accounts for the 
correlation of the electron at fee = bound with holes 
in other kh states. Excitons with all values of the cen- 
ter of mass wavevectors are involved in this peakja as we 
will also show in the next section. This correlation struc- 
ture is of course not present in the Born approximations, 
which maintain their single quasi-particle peak structure 
at any density. As the density rises, the relative weight 
of the satellite structure in the SCLA spectral function 
with respect to the main quasi-particle peak increases. 
Moreover, both structures become broader. For densi- 
ties n > 0.1, it becomes diflficult to distinguish between 
the excitonic and the main peak, and a single broad, red- 
shifted quasi-particle structure appears. The merging of 
the excitonic peak with the single quasi-particle peak can 
be interpreted as a Mott transition. We remark how- 
ever that it is difficult to give a rigorous estimate of the 
Mott density, because broadening is gradual and a sharp 
threshold is not evident. 
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FIG. 4. Self energy E(fc = 0,lu) for T = 0.1 and differ- 
ent densities given in the Figure, in the self consistent Born 
approximation (SCBA) and self-consistent ladder approxima- 
tion (SCLA). 

In Fig. y we show the imaginary part of the electron 
self energy for different densities. We will also refer in the 
following to this quantity as (energy dependent) broad- 
ening. In the SCLA 3(S) shows two peaks corresponding 



to dephasing experienced by an unbound or a bound state 
propagated in the system. Both peaks increase with den- 
sity. At low density, n = 0.01, dephasing at the exciton 
energy dominates, while at n = 0.05 both peaks become 
comparable. In the SCBA, the peak of 3(S) at — -Ef, is 
obviously absent, as propagation of bound pairs is not 
allowed in the theory. Instead at low density n — 0.01, 
the broadening at the main quasi-particle energy {lo — 0) 
is similar in the SCLA and SCBA. Indeed, few excitons 
are expected in the plasma, and scattering mainly origi- 
nates from free carriers. An heuristic understanding may 
be obtained with the action mass law relating the con- 
centration of different chemical species in a reaction. In 
our case, electron-|-hole ^-> exciton, and the law states 
nc^/nx — n*(T), where ric and nx represent the den- 
sity of unbound carriers and excitons respectively, n* (T) 
is a cross-over density which depends only on the tem- 
perature, and at T=0.1, n* ~ 0.005 . For n < n*, 
nx <C n. Therefore, a small density of excitons is ex- 
pected a.t n = 0.01 ~ n*. At this density the exciton-free 
carrier contribution to free carrier broadening is weaker 
than the free carrier-free carrier contribution. At larger 
density, bound excitons in the plasma become dominant, 
and broadening at the quasi-particle peak becomes much 
larger in the SCLA than in SCBA. This trend is observed 
up to n = 0.15 indicating relevance of correlated states 
in the plasma even at this large density. However, the 
heuristic interpretation of the broadening based on the 
action mass law given above becomes meaningless, as it is 
not possible to distinguish between two chemical species 
anymore. At the highest considered density n — 0.4, we 
clearly observe a dip to very small broadening in the self- 
energy for SCBA case. This dip is exactly at the Fermi 
energy and accounts for the blocking of the scattering 
at the Fermi level. In fact, only at this large density the 
Fermi gas becomes degenerate (Fermi energy much larger 
than T). In the SCLA, the broadening never vanishes, but 
is only partially reduced at the Fermi level. This is in- 
dicative of a more complicated structure of the electron- 
hole plasma and of the broadening process, where Pauli 
blocking looses much of its effectiveness, and is suggestive 
of a non- Fermi liquid behavior. We finally noticed that 
3?(E) is comparable to S(S) at any density in the SCLA, 
and thus remarked that the excitonic satellite peak in 
the spectral function does not correspond to a zero of 
w — 3fi(S(u;)). Indeed, we do not expect appearance of an- 
other simple quasi-particle in the plasma. For w ~ —Eb, 
and for Et » |3(E)|, we have Ak^o -- \'^{J:)\/\Eb\^. 
Thus, the correlation satellite in the spectral function 
follows the peak of 3(S) at uj = — £"&. 

In Fig.^ we show the density dependence of the elec- 
tron (or hole) chemical potential as a function of the den- 
sity for T = 0.1, in the three considered approximations. 
At very low density, n <C 0.01, the chemical potential 
is fi <ti —0.4, and its value is similar in all approxima- 
tions. In this case, we are describing free electron-hole 
pairs, as also suggested by a law of mass action, which 
gives a crossover density of about n* — 0.005 at T=0.1. 



In this case the description of /i is rather good even in 
approximations that neglect existence of bound states. 



0.0 



-0.1 






o 

D- 



-0.2 



E 

0) 

O 



-0.3 - 



-0.4 



-0.5 











^ 










-^'^ '^J'^ 








^^'^'^ 


,^^,X'^ 


- 




,^^--'\^ 


^ 


- 


/ 

/ 

/ 


;; 


'e^-- 










i// 

1 




SCLA 

Bom 

SCBA 








1 


bosons 

1 


1 





0.1 



0.3 



0.4 



0.2 

Density 
FIG. 5. Chemical potential ^ in the three approximations, 
at T=0.1. The short-dashed line shows half the chemical po- 
tential of non-interacting bosons, with ground state energy of 

However, for n > 0.01, the chemical potential in the 
Born approximations is much larger than for the SCLA. 
In fact, in the SCLA we are also describing the fraction 
of cold interacting excitons in the plasma. As the ground 
state energy of this part of the plasma is at about —Et, 
the total energy of the system is thus reduced with re- 
spect to that of a gas of unbound particles. In order to 
clarify further this point, we plot on the same graph half 
of the chemical potential of bosons at a ground state en- 
ergy of —E}, and mass equal to the exciton mass. This is 
the chemical potential of a gas of electron (and holes) 
completely bound into bosonic excitons. The halving 
comes from the equilibrium condition ^x = Me+M/i = 2/.t. 
The chemical potential of bosons compares reasonably 
well with the chemical potential calculated in the SCLA, 
up to /Lt < —0.25, n < 0.1. Above this limit, the chem- 
ical potential from the SCLA grows faster. There are 
two reasons for this faster growth: first, the exciton gas 
is repulsively interacting, second, the exciton gas even- 
tually undergoes a Mott transition and ionizes into free 
carriers. At densities n > 0.2, the SCLA and SCBA are 
indeed comparable, and at even larger density, also the 
MBA approximation is reasonable. This shows that the 
SCLA is a powerful tool for the investigation of the inter- 
mediate regime of densities, where deviations from both 
the (non-interacting) bosonic and pure fermionic models 
arc important. 



IV. MAPPING THE SCLA TO A BOSONIC 

MODEL AT LOW DENSITY AND 

TEMPERATURE 

In Fig§ we plot A{k,uj) for T=0.043 and n=0.02. We 
can observe the parabolic dispersion of the main quasi- 
particle peak, broadened by the scattering (white region 
in the plot). Broadening is larger at small k due to the 
larger phase space in one dimension. For w < 0, we have 
the correlation structure shown for fc = in Fig. pi which 
is also present at larger k. We remark that this region 
of correlated electrons extends into a region of k of the 
order of 1 (i.e. of a^^). Its shape (dispersion) is related 
to the exciton wavefunction, as we show in the following. 
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FIG. 6. Contour plot of Ae{k,uj) at T = 0.04 and n=0.02. 

In a low-density, low-temperature limit, we can start to 
consider G = Gq, and analytically calculate the retarded 
pair propagator H^{k,Lu) from Eq. (p^, neglecting Gq , 
which is of the order of the density. We obtain: 
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where 7 > is the usual regularization number. S olvin g 
for the T-matrix with the Bethe-Salpeter equation (14a), 
and expanding around its pole we obtain: 
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Thus, the T-matrix at the lowest order for the contact 
potential has the form of a free propagator, for particles 
at an energy —1/2 + fc^/8, having a mass which is twice 
the electron mass, i.e. the electron plus hole mass. Ex- 
panding the Bethe-Salpeter equation to higher order, we 
obtain: 

r(i)+(fc,w) = T(°)+(fc,tj)ij(i)+(fc,cj)T(°)+(fc,tj). 



Therefore, we may interpret T^^^~^{k,uj) as the free exci- 
ton (boson) Green function, and H^^'>^{k,uj) as the low- 
est order self-energy.Ell In particular 



ij(i)(l,2) = 2iG(i)(l,2)Go(l,2), 



(16) 



G'^^\k,oj) = Go{k,uj)J:^"\k,uj)Go{k,uj), 

S]("ni,2) = iT(")(l,2)Go(2,l). (17) 

As T^°)(fc,w) is peaked close to w = ^1/2, we may ne- 
glect GQ{k,Lu) in this region, and use 
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Integrating around lo' = —1/2, we obtain: 
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where n is the total density. This is the structure shown 
in Fig. U at a; = —1/2, which produces the satellite 
correlation peak in the electron spectral function, shown 
in Fig. pi We also notice that this structure has negative 
dispersion, as also apparent in Fig. |6|. 
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FIG. 7. The perturbation expansion to lowest order for the 
boson (exciton) self-energy, and the four-particle vertex _F'°' . 
Thin lines are the bare electron (or hole) propagator. The 
shaded two-particle Green function is the bare exciton prop- 
agator. 

We can now define an exciton-exciton interaction at 
low density, from the Eqs. (Qq-O) above, which are pic- 
torially shown in Fig. |7|. We thus define a four-point in- 
teraction kernel i^'^"^ and the boson self-energy iJ'^^^( 1,2) 
is written as: 

ij(i)(l,2)= /dl d2F(")(l,2;l,2)T(°)(2,l). (18) 



We can see from Fig. R that i^^"' represents the electron 
(or hole exchange) in the scattering of the two excitons. 
Thus 

i^(")(l, 2; 1'2') = 2iGo(l, 2)Go(l', 2')Go(l, 2')Go(l', 2). 

(19) 



If we neglect retardation (or memory) effects, we may 
define an instantaneous potential at the exciton energy 
using Ti = T2,Ti> — T21, and calculating the resulting F at 
the exciton frequency wq = —1/2 — Eb- For simplicity, 
we also use k,k' '^ for the incoming excitons, and arrive 
to 
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Wo J 



2z^G, 
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-w')Go+(g'-9,7--'') = 
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(20) 



642V4 + g2-(4-g2) 



2\2 



92 (4 + ^2) 

This function rapidly drops to zero at g > 1 as expected, 
and F^^\0,ujq) = 6. This value is also obtained from 
standard boson-boson exchange-interaction expressions 
(see e.g. Ref. ^2|) when only the electron-hole interac- 
tion is used. We show the full structure of F^-^^q^uJo) In 
Fig. g 
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The exciton-exciton interaction potential 
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We are now in a position to discuss some of the qual- 
itative changes introduced by simplifying the Coulomb 
interaction to an on-site one. In the limit of low den- 
sity, we may in fact compare the exciton-exciton inter- 
action calculated above with that calculated using a re- 
alistic Coulomb interaction, as a standard boson-boson 
exchange expression at small momenta existsH3 In a re- 
alistic wire, a short-range cut-off is introduced by the 
finite size of the electron and hole wavefunctions in the 
confinement directions. A typical cut-off is of the order 
of the Bohr radius or smaller. Typically, we have to con- 
sider tight confinements in order to avoid participation of 
higher confined levels into the exciton wavefunction. The 
Coulomb interaction is then reasonably represented by a 
function Vr = A/ [r -I- ro), where ro is a cut-off distance, 
and A is then chosen such that Ef, = —1/2. We used 



To = 0.1, and A = 0.18. In this case the resulting Bohr 
radius is 1, indicating tight confinement of the carriers 
in the confinement plane. The exciton-exciton exchange 
interaction at small q is calculated in the small q <^ as^^ 
limit 



W^^ 



2^Vr{k 

k,k' 



k')(t,u{k)c^is{k') 



[\cbis{k)\'' - q^is{k)^is{k')] -1.4 



(21) 



where (j)is{k') is the Is wave- function of the exciton. The 
positive term in Eq. ^ is the exchange term due to 
the attractive electron-hole interaction, while the neg- 
ative term is due to the electron-electron and hole-hole 
repulsive interactions. In the delta like-potential the neg- 
ative term is absent because of locality. In the long range 
case instead, the electron-electron (and hole-hole) inter- 
action largely cancel the electron-hole interaction in the 
exchange integral, resulting in an exciton-exciton inter- 
action which is four times smaller when compared with 
the value of 6 obtained for F^'^^q = 0,a;o). We conclude 
that local interaction leads to an over-estimation of the 
boson-boson interaction at a given density, and therefore 
of broadening. However, this does not imply that broad- 
ening effects are qualitatively different. Moreover, we also 
stress that this overestimation does not concern at all the 
comparison of the SCLA with the other simpler approx- 
imations considered in this work, as all are carried out 
using the same interaction potential. 



V. OPTICAL PROPERTIES 

We now consider interaction of the electron-hole sys- 
tem with a transverse electromagnetic field in the dipole 
approximation: 
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where 



Cq.q — 



I At: he 
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rriocV qtV 



^11 (q), 



is the dipole matrix element, which includes also the over- 
lap integral of the electromagnetic field with the con- 
fined carriers /y . Here q represents 2D wavevectors in 
the plane orthogonal to the wire axis, q are wavevectors 
along the wire, while qf = q^ + q"^, mo is the free elec- 
tron mass, c the velocity of light, and Pcv the momentum 
matrix element between the conduction and hole bands, 
flq^q is the photon destruction operator, and Pq is the po- 
larization operator, defined as Pq = ^j, dq-kCk, i-e. local 
in real space in the dipole approximation. As the section 
of the wires is always much smaller than the wavelength 
of light, we have /|| — 1. Moreover qt — Eqap /{he) is 



almost constant, thus we use arbitrary units from now 
on, with C = 1. Actual values of absorption and gain are 
easily calculated by using the appropriate value of C for 
the given material. 

In the process of scattering of light from the system, 
an absorbed photon creates a coherent electron-hole pair 
which propagates in the crystal, interacts with the back- 
ground plasma, then recombines, and a photon is emit- 
ted. When we introduce photon propagators, the interac- 
tion process is represented by a photon self-energy. We 
will not solve the photon Dyson equation, nor we will 
dress the electron and hole propagation with the photon 
interaction. In other words, we neglect both polaritonic 
effects, and higher order non-linear interactions. In fact, 
we are restricting our attention to weak external fields, 
i.e. assuming interaction of the plasma with the photon 
field to be weak. Otherwise, it could be possible for the 
system to be driven out of equilibrium, because the pho- 
ton chemical potential is zero and far below that of the 
elecron-hole plasma. This equilibrium condition is well 
fulfilled in real systems, exception made for the strongest 
laser fields. Within this weak interaction approximation, 
the photon self-energy is exactly given by the pair corre- 
lation function 



n,(r,r') = (TPt(^)P^(^')). 



(23) 



However, this correlation function has to be calculated 
within some approximation, and also using approximate 
propagators. The resulting self-energy may not be con- 
serving. In particular, particle flux conservation at the 
vertex translates in the longitudinal f-sum rulellj, which 
has to be respected in any reasonable approximation. 

In order to shed light on this important issue, we 
consider the problem of interaction of photons with the 
plasma from a different viewpoint. As we are only inter- 
ested in the linear response to the external electromag- 
netic field, the correct photon self-energy is also given by 
the linear response of the plasma to a classical electro- 
magnetic field, i.e. considering the photon creation and 
destruction operators in Eq. ( p^ ) as c-numbers. Baym 
and Kadanoff give a recipe to construct a conserving ex- 
pression for this response function in Ref. ttq. When 
an external coherent electromagnetic field is applied to 
the plasma, coherence develops in the plasma, and pair 
propagators have additionally to be defined along with 
the single particle propagators. They are analogous to 
anomalous propagators in the theory of superconductiv- 
ity. The electron-hole anomalous propagator reads 



Ge,.(l,2) = z(rd(2)c(l)). 



(24) 



Times are defined on the Keldysh contour as usual. An- 
other anomalous propagator Ghe can be similarly de- 
fined for convenience, and notation regrouped into a ma- 
trix one, like for Mambu propagators in the problem of 
superconductivity.cj Dyson Eq. (0) is extended to in- 
clude both anomalous propagators, and the external in- 



teraction H„ 



d2 G(7i(l,2)G(2,3)= (5(1,3)+ / d4 E(1,4)G(4,3) 



+ / d4 Eeh(l, 4)G,,e(4, 3) - A{l)Ghe{l, 3) (25) 



d2 G^\l,2)G,h{%i) ^ / d4S(l,4)Ge/,(4,3) 



+ / dl Se/,.(1, 4)G(4, 3) ~ A(1)G(3, 1) 



(26) 



As a general rule for generating a conserving approxi- 
mation, the self energies have to be obtained from func- 
tional derivation of a grand-potential with respect to the 
propagators. The normal self-energies have been already 
introduced. For the anomalous self energy Se/i(l,2), we 
need to generalize the expression of the grand-potential 
including terms where the normal propagators are re- 
placed by the anomalous ones. The anomalous propaga- 
tors Geh are at least first order in the external potential 
A, as spontaneous symmetry breaking is ruled out by 
the Mermin- Wagner theorem. For linear response, we 
need the anomalous self energy Yi^h up to linear terms 
in the external potential only, thus at most linear terms 
in the anomalous propagators. Consequently, we need at 
most quadratic terms in these propagators in the grand- 
potential to obtain the Sg/i from the functional deriva- 
tive. We conclude that the only and sufficient term to 
be considered is the anomalous Fock term, which is the 
usual Fock term with anomalous propagators replacing 
normal ones. As the approximation is conserving, we as^ 
also guaranteed that absorption fulfills the f-sum ruleEJ 
for any value of density and temperature. The anoma- 
lous self-energy obtained from the anomalous Fock term 
reads 



Seh(l,2) = -iGe,.(l,2)-HO(A2). 



(27) 



Normal G are at least second order in the external po- 
tential A, while both Geh and T,f.h are at least first order. 
The anomalous Dyson Eq. (Efl) to lowest order becomes: 



d2 G-i(l,2)Ge/.(2,3) = -i dl Ge„(l,4)G(4,3) 



-A(1)G(1,3). 
Introducing the response function 

6G,h{l,l) 



n(l,2) = -* 



A{2) 



(28) 



(29) 



A=0 



and taking the first order variation of Eq. (Eq), we obtain 
and equation for 11(1, 2): 

11(1,2) =iG(l,2)G(l,2)-i /d3G(l,3)G(l,3)n(3,2) = 



(30) 



= i?(l;2)- f d3 H{l;3)n{3,2), 
which is easily solved and gives 
n(l,2) =iJ(l;2)+ / d3d4iJ(l;3)T(3;4)iJ(4;2). (31 



It is straightforward to show that 11(1,2) obeys bosonic 
Kubo-Martin-Schwinger relations, as both H and T do. 
In particular, we have: 



n+{q,Lu)=H+iq,Lu)\T+iq,u)\ 
Il<{q,Lo)^H<iq,Lo)\T+iq,Lu)[ 



(32) 
(33) 




Absorption a(w) and photoluminescence PL{uj) at q 
are derived as a{uj) ~ — 3[n+(0,a;)] and PL{llj) = 
— 9[n<(0,w)] respectively. 

We remark that higher order terms than the anomalous 
Fock one have to be considered in the anomalous self- 
energy expansion when a finite external field is present, 
such as in the case of a strong laser field. Anomalous 
Born terms calculated in the Markov approximation have 
been considered in the semiconductor Bloch equations by 
Lindberg and Koch in Ref. Ea, and named polarization- 
polarization scattering terms. Here we do not further 
pursue such extensions, which are clearly beyond the 
scope of the present paper. We only remark that our ap- 
proach is easily extended to such conditions, while keep- 
ing full control of the conservation laws. 



A. Optical spectra and stability of the excitonic 
emission 

We show in Fig. |[ the absorption spectra in the nor- 
mal direction {q = 0) at T = 0.2, for different densities, 
in the MBA, SCBA, and SCLA. In the absorption spec- 
trum, we remark at low density the characteristic exci- 
tonic peak and the continuum absorption. The exciton 
linewidth is extremely narrow in the SCBA, and much 
broader in the MBA. The SCLA linewidth is interme- 
diate. In the SCBA there is a very small broadening 
at the exciton energy, as no excitons are represented in 
the theory, while in the MBA, the free carrier broad- 
ening at u; = is assumed. This is somewhat larger 
than the broadening in the SCLA at T = 0.2. For large 
density n > 0.2, the excitonic peak is bleached in all ap- 
proximations, and a region of negative absorption, i.e. 
optical gain, appears. The absorption changes sign at 
w = /ie -I- /ift = 2/i, coinciding with the change of sign 
of the Bose factor, and resulting in a well defined -i.e. 
positive- photoluminescence through the Kubo-Martin- 
Schwinger relation. 



We plot the emission or PL spectra in Fig. 10, First we 
notice the different scale for the MBA, where a smaller 
peak emission is calculated even for the largest consid- 
ered density. For large density n > 0.2, we observe a 



saturation of the intensity. This can be understood in 
terms of the fermionic nature of the carriers, when the 
electron- hole plasma becomes degenerate. Second, we 
notice exciton-like emission even when excitonic absorp- 
tion is bleached in all models. In fact, vertex corrections 
place the free carrier emission at the exciton energy, even 
when no bound excitons are described in the gas. 
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FIG. 9. Absorption spectra in the Markof-Born Approxi- 
mation (MBA), self consistent Born approximation (SCBA) 
and self consistent ladder approimation (SCLA) for different 
densities indicated in the figure. 

Eventually, only at very small density the free car- 
rier and excitonic emission become distinguishable, as we 
notice an exponential emission shoulder at high energy 
w > in the SCLA and SCBA, reminding of a fermionic 
emission tail at small degeneracy (Boltzmann distribu- 
tion). At larger density, the two features at different 
energy merge into a unique one, due to increased broad- 
ening. A single peak is also observed in experiments. 
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where finite noise in the data and Large inhomogeneous 
broadening mask any minor feature.cZI For this reason, it 
is in practice very dehcate to estabhsh the position of the 
band-gap (bottopi of the free-carrier bands) directly from 
emission data.E3 Third, we remark that the low-energy 
shoulder of the PL emission in the MBA is Lorentzian 
because of the Markov approximation, and the broad- 
ening largely over-estimated also at low energies. The 
linewidth in the SCBA is instead underestimated at low- 
density, as for absorption. 




spectral function. Interestingly, exciton bleaching in ab- 
sorption appears at a density which is comparable to that 
where the band gap renormalization crosses the emission 
energy, i.e. n ~ 0.15. This can therefore be assigned as 
a Mott density, i.e. the energy where the binding energy 
becomes negligible with respect to broadening, thus ef- 
fectively vanishing. There are two interesting features to 
be observed in Fig. 11: first, the band gap renormal- 
ization is negligible in the SCLA for n < 0.1, second, 
emission energy is constant in this range of densities, 
and blue-shifting less than in the SCBA for larger den- 
sity. Comparison with the MBA is vitiated by excessive 
broadening in this approximation. The simple Hartree- 
Fock approximation instead gives results which are sim- 
ilar to those obtained in the SCBA for the band gap 
renormalization and the PL jernission peak, apart from 
more pronounced blue-shifts.E£l Stability of the emission 
peak is usually interpreted as a partial compensation of 
the self-energy and vertex corrections. As in the Hartrce- 
Fock approximation the broadening effects are missing, 
we deduce that these effects are indeed relevant for the 
cancellation found in the SCBA and in the SCLA at high 
density, and explain the reduction of the blue-shift with 
respect to simple Hartree-Fock calculations. 




FIG. 11. The band gap renormalization (BGR) and emis- 
sion peak energy in the SCBA and SCLA as indicated in the 
figure. 



B. Excitonic gain 



FIG. 10. Photoluminescence spectra in the MBA, SCBA 
and SCLA for diff'erent densities indicated in the figure. 

An interesting -and directly observed- physical quan- 
tity, is the energy of the PL peak. As remarked above, 
emission is excitonic due to the correlation of the coher- 
ent electron-hole pair emitted. In Fig. |ll| we plot the 
position of this peak in the PL as a function of the car- 
rier density, together with the band-gap renormalization 
(BGR) defined as twice the energy of the main peak of the 



Optical gain is usually considered in the regime of de- 
generate electron-hole plasma. In this regime, spectral 
functions can be assumed to be simple Lorentzians, and 
the absorption can be written as: 



a{uj) ex 



dfc 1 - 2/(^2/4 - /z) 



2tt (w- 



fc2/2)2+^2^ 



(34) 



where 7 is the spectral broadening. Then, negative ab- 
sorption or gain occurs only when the Pauli blocking fac- 
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tor (1 — 2/(fc^/4 — n)) becomes negative, which necessar- 
ily requires a chemical potential above the band gap, or 
in other words, inversion. However, both vertex correc- 
tions, and deviations of the spectral function from simple 
Lorentzian should be taken into account at lower density 
and temperature. In this case, absorption can be written 
as: 



, , ,^^ , , ,^ f dk du! 



(35) 



[1 - f{^ -^'- /^) - f{^' - A^)] Akiu: - ij')Ak{iu') , 

The term in square brackets 1 — /(cj— a;' — /i) — /(cj' — /i) is 
the generalization of the Pauli blocking factor for a sys- 
tem where the quasiparticles are described by arbitrary 
spectral functions, and clearly one of the Fermi functions 
refers to electron occupation, the other to hole occupa- 
tion. This expression of the absorption is valid only in 
the case of short range potential, and for the long range 
case T"*" depends also on the relative momentum of the 
pair k, and becomes part of the integral kernel. Using 
the definition of Fermi and Bosc functions, the absorp- 
tion can be recast in the form 



a{uj) 






dk duj' 
27r^ 



(36) 



Xf{uj - w' - /i)/(^' " fJ')Aki^ - Uj')Ak{Lj') , 

where 17(0; — 2/i) is the Bose function. Therefore, gain 
clearly occurs below twice the chemical potential /i (i.e. 
fie + fJ'h), even when the chemical potential is below the 
band-gap. However, for it to be sizeable, the spectral 
functions must have non-negligible weight in this region 
of the spectrum. In the MBA and SCBA case, this weight 
is clearly given by broadening effects alone, while in the 
SCLA it is also due to the presence of the excitonic corre- 
lation peak in the spectral function as shown in Fig.pl A 
sizeable enhancement of the gain is found in correspon- 
dence to the excitonic resonance because of the vertex 
correction, given by the \T^\'^ factor in Eq. (Bq). When 
the 2fi < —Eb, both gain and excitonic absorption coex- 
ist. In this case, we talk of excitonicgain. The system is 
inversion-less in the sense of Eq. ( p4| ) above. However, 
only in the SCLA we are describing gain due to the pres- 
ence and scattering of excitons in the low temperature, 
low-density plasma, while in the other approximations, 
gain is purely related to dynamical effects in the interac- 
tion vertex with the photon. 

In Fig.|l^ we plot the absorption spectra close to the 
exciton resonance at T = 0.1 <C Eb- For the considered 
densities n < 0.1, gain coexists with the excitonic reso- 
nance in absorption. Only for the MBA at n = 0.1, the 
absorption peak is completely shifted to higher energies, 
and we can not anymore speak of excitonic absorption. 
As usual, the broadening in the exciton spectral region 



is too small and unphysical in the SCBA, resulting in 
very sharp and unphysical features. Gain in the SCLA 
is instead over a larger spectral region, and its value is 
larger than what predicted in the MBA at large density 
n > 0.3 shown in Fig.0. Excitonic gain clearly originates 
from inclusion of excitons in the plasma for the SCLA. 
In fact, reasonable broadening is calculated in the exci- 
ton spectral region, and mainly originating from exciton- 
exciton scattering as shown in Sec. IV. We conclude. 



that the SCLA indicates that sizeable excitonic gain can 
be obtained at moderate density and temperature, in a 
reasonably large spectral region of the order of a fraction 
of Eb. 



n=0.03 
n=0.05 




FIG. 12. Absorption spectra at T = 0.1 for the MBA, 
SCBA, and SCLA, and low densities indicated in the figure. 

Excitonic gain has been widely studied in II- VI quan- 
tum wellsEZI. However these are two-dimensional system, 
and excitonic correlations are expected to be less pro- 
nounced in this systems, so at this stage it is not reason- 
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able to make even a qualitative comparison. For quantum 
wires, sizeable gain at 10 K (T ~ 0.05) for an estimated 
carrier density well below the Mott has been recently 
claimed by Sirigu et alE3 This is a first indication that 
excitonic gain might be relevant in these systems. How- 
ever, a qualitative and quantitative understanding of this 
experiment is clearly premature. From the experimental 
side, inhomogeneous broadening due to interface disorder 
has to be decreased well below the binding energy, while 
a more realistic Coulomb interaction has to be addressed 
by the theory. 



VI. CONCLUSIONS 

We have presented a model that includes excitonic 
correlation in the description of a highly-excited semi- 
conductor quantum wire. The model has been sim- 
plified using a short range potential, and considering 
a polarized gas. Correlation has been calculated self- 
consistently at the ladder level (SCLA). We have shown 
that bound states appear as a low-energy correlation 
peak in the spectral function of electrons and holes. We 
compared the results obtained with SCLA to the ones ob- 
tained within lower order approximations that do not in- 
clude excitonic correlation (Born approximations). Even 
though the model is purely fermionic, we have shown 
how it can be effectively mapped at low temperatures 
and density to a gas of interacting excitons. We have 
thus derived an analytical expression for the effective 
exciton-exciton interaction. The linear optical properties 
of the system have been calculated including vertex cor- 
rections at the Fock level. This ensures the conservation 
of sum rules in the optical response. The excitonic ab- 
sorption at low density ensues both in ladder and Born 
approaches, but the broadenings at the exciton energy 
are either too small or too large in the Born approx- 
imations, depending on whether frequency dependence 
of the broadening is included or not. Excitonic emis- 
sion well beyond exciton bleaching is also predicted in 
all models, but the peak shifts are more pronounced in 
the Born approximations, showing that a better cancella- 
tion between the self-energy and vertex correction results 
from the introduction of excitons. We have also shown 
that sizeable excitonic gain can be predicted at low tem- 
perature and density, when the electron plus hole chem- 
ical potential is just below the exciton energy. However, 
its correct description must include exciton broadening 
from exciton-exciton scattering, and unphysical values 
are thus obtained in the Born approximations. These 
qualitative conclusions clearly hold even for more refined 
descriptions of the electron-hole plasma than the Born 
approximations considered here, when these descriptions 
do not account for excitons in the plasma. For example 
this is the case for the model recently presented by Das 
Sarma and Wang in Ref. ^, which includes screening at 
the plasmon-pole approximation level, but no excitons in 



the electron-hole gas. The excitonic gain calculated us- 
ing this approximation shows a divergent behavior close 
to —Eb, similar to what we have found for the SCBA. 
The extension of the proposed approach to a full non- 
equilbrium theory could contribute to a deeper under- 
standing of many body effects in electron-hole systems, 
also in the perspective of ultrafast optical experiments 
than can investigate far- from equilibrium conditions. In 
this context a recent theory by Hannewald et alEi pro- 
vides a convincing step in this direction, yet remaining 
within a non-dynamical approach. The last issue concern 
the intrinsic limitations of the SCLA. It has been shown 
that in order to describe correctly the very diluite exci- 
tonic limit additional diagrams must be added to the self 
consistent ladder approach. An analysis of this problem 
can be found in Ref. pll where a generalized T matrix 
approximation for the fermionic self energy overcoming 
this problem has ben proposed. In conclusion, we have 
clearly shown that in the definite and important physi- 
cal region oi T < E^ inclusion of excitonic correlation in 
the electron-hole plasma is relevant and necessary, and 
that the simplified model presented in this paper is a 
well-understood starting point for this purpose. Thus, 
its further developments to address realistic systems in 
higher dimensions are well motivated. 
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